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Abstract: As a step towards understanding the properties of string theory in time- 
dependent and singular spacetimes, we study the divergence of density operators for string 
ensembles in singular scale-invariant plane waves, i.e. those plane waves that arise as the 
Penrose limits of generic power-law spacetime singularities. We show that the scale invari- 
ance implies that the Hagedorn behaviour of bosonic and supersymmetric strings in these 
backgrounds, even with the inclusion of RR or NS fields, is the same as that of strings in 
flat space. This is in marked contrast to the behaviour of strings in the BFHP plane wave 
which exhibit quantitatively and qualitatively different thermodynamic properties. 
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1. Introduction 

The study of string theory in the cosmological setting of non-trivial time-dependent and 
possibly singular backgrounds is an outstanding problem. While string theory in general 
time-dependent backgrounds is hard, in principle string theory in plane waves (or, say, flat 
orbifolds) is exactly solvable. Hence one can hope to gain some insight into the nature and 
properties of string theory in general time-dependent backgrounds by studying it in the 
context of time-dependent plane waves (or time-dependent orbifolds). With this motiva- 
tion, in Ijl], |2[ ^, ^ we set out to systematically identify tractable yet physically interesting 
time-dependent plane wave metrics and develop tools for the quantisation of string theory 
in such backgrounds.^ 

In m, 1^ it was found that the plane wave metric 

ds^ = 2dudv + Aab{u)x''x^du^ + dx^ (1.1) 

associated to every space-time metric and choice of null geodesic in that space-time by 
means of the Penrose limit construction ^, encodes covariant information about the 
rate of growth of curvature and geodesic deviation along the null geodesic. In particular, 
singularities of Aab{u) result from curvature singularities of the original space-time. 

Moreover, in [^, Q we studied the Penrose limits of metrics with space-time singu- 
larities. We found that the resulting plane wave metrics for a large class of black hole, 

^For a detailed review of recent work on time-dependent orbifolds see [0]. 
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cosmological and null singularities (namely all metrics with singularities of power-law type 
g, |10|, |ll| subject to some energy condition) exhibit a remarkably universal leading 



behaviour, 

Aab{u)^u-'' , (1.2) 

near the singularity. Plane waves with such a wave profile are scale invariant, i.e. invariant 
under the scaling (boost) (n, v) —>■ (An, A"^^) and hence also homogeneous j^, |l2|, This 
scale invariance will turn out to play an important role in the following. 

Thus the considerations of ^] single out the singular scale-invariant plane waves with 
profile ~ u~'^ as the backgrounds to consider in order to obtain insight into the properties 
of string theory near physically reasonable space-time singularities. String theory in these 
singular homogeneous plane wave backgrounds is exactly solvable 1 13 , |l^, Q , and various 



aspects of string theory in this class of backgrounds have already been studied in particular 
in in 



Recently, motivated by the BMN correspondence |14], thermodynamical aspects of 
string theory in the BFHP plane wave |l^, |l^] and other plane waves with constant 

(i.e. u-independent) profile Aab, have been investigated in e.g. |l^-[^. In particular, it 
was found that the Hagedorn behaviour of strings in such backgrounds (and with RR flux) 
differs quantitatively and qualitatively from that of strings in flat space. 

Here we wish to study the analogous question for strings in the singular scale-invariant 



plane wave backgrounds with profile (1^). In order to address this issue we need to first 
come to terms with the fact that string theory in these time-dependent singular homoge- 
neous plane waves leads to a time-dependent light-cone Hamiltonian. Thus the study of 
the "thermodynamics" of such systems requires some care. 



Following |2g, |27|, one can study the evolution of non-equilibrium systems with 
respect to invariants of the system. Given such an invariant /, i.e. an operator which 
satisfies dl/dt = 0, one can introduce an analogue of the Boltzmann thermal state, the 
density operator 

PI = Jot ■ (1-3) 

This density operator satisfies the quantum counterpart of the classical Liouville theorem 
for the phase space density, namely tr/5/ = 1 and dpi/dt = 0. A convenient choice of 
invariant for a system with a time-dependent Hamiltonian H(t) is I = H(to), as it reduces 
to the standard choice in the case of a time-independent system and provides an "adiabatic" 
approximation to the system provided that H{t) varies sufficiently slowly with time near 

t = to. 

We will study the corresponding "thermal" partition function 

Zj{(3) = tre-f^^ (1.4) 

for strings in the singular scale-invariant plane waves. Our main result is that the partition 
fuction for bosonic or type II strings, even with the inclusion of RR fields, diverges at a 
critical value Pc of the parameter /3 which is identical to the inverse Hagedorn temperature 
Pc = l/Tf/ of strings in fiat space |2^. This is in marked contrast to the thermodynamical 
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behaviour of strings in the non-singular constant Aab plane waves which, as mentioned 
above, can be quite different. 

The calculation establishing this result highlights the significance and implications of 
the scale invariance of these plane waves. Indeed scale invariance can be seen to imply 
that the string mode frequencies are independent of the light-cone momentum and hence 
uniformly approach the flat space frequencies at large string mode number n. 

This may be an indication that string propagation in these scale invariant plane waves, 
which arise as Penrose limits of space-time singularities, has properties more in common 
with the propagation of strings in flat space than either have with strings in the BFHP 
and other time- independent plane waves. 

In section 2 we summarise the approach of p6| , |2^, to non-equilibrium thermody- 
namics via invariants. In section 3, following [^] we present the light-cone Hamiltonian of 
type IIB superstring theory in singular homogeneous plane wave backgrounds and exhibit 
its relation to our preferred invariant. The implications of the scale invariance are explored 
in section 4 and the detailed computation of the partition function is given in section 5. 



2. Invariants and Thermodynamics of Time-Dependent Systems 



For present purposes we will adopt the point of view (see e.g. [gg, |28|]) that suitable 
analogues of the Boltzmann thermal state for a time-independent Hamiltonian system, 

PH = loTT , (2-1) 

can be constructed as density operators of the form 

Pi = jTf , (2.2) 

tre-/'^ 

where / is an invariant of the system, i.e. a possibly explicitly time-dependent operator in 
the Heisenberg picture satisfying the equation 

1/= + ,1=0. (2.3) 

Note that for a time-independent system evidently the Hamiltonian H itself satisfies this 
equation. In particular, the density operator satisfies the quantum counterpart of the 
classical Liouville theorem for the phase space density, namely tr /5/ = 1 and the Liouville- 
von Neumann equation 

J,Pi = §IPI + ^iH,M=0 . (2.4) 

As a consequence of this equation, the density operator pi{t) allows one to calculate the 
time evolution of the expectation value of any operator O in the mixed state characterised 
by the invariant I, 

<0>i (t) =tv{Opi) . (2.5) 
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Different cfioices of / correspond to different initial "thermal" ensembles with partition 
functions 

Ziip) =tre~^^ . (2.6) 



In particular, in |^6|] it has been shown that this reproduces and unifies various approaches 
to non-equlibirum thermodynamics such as mean-field or Hartree-Fock methods. 

The parameter f3 is in many ways analogous to an inverse temperature 1/T. However, 
here and in the following we occasionally put the word "thermal" in quotes to emphasise 
that we are not claiming that (2.2) describes "the system at temperature P = 

For the harmonic oscillator systems under consideration (namely the light-cone Hamil- 
tonians for strings on time-dependent plane waves), invariants are easy to come by, e.g. 
using invariant oscillators (Appendix A).^ 

A convenient, but by no means the unique acceptable, choice is the "instantaneous" 
density operator 

/ = H{to) (2.7) 

(which will be explicitly time-dependent when written in terms of Heisenberg operators). 
The corresponding density operator reduces to the standard choice in the case of a time- 
independent system, and it provides an "adiabatic" approximation to the system provided 
that H{t) varies sufficiently slowly with time near t = tQ. 



3. The Light-Cone Hamiltonian and the Invariant 

We now summarise some results from [jl2| regarding the light-cone Hamiltonian for bosonic 
strings in a purely dilatonic singular homogeneous plane wave background (we will briefly 
discuss other string theories and backgrounds at the end of section 5). 
In the light cone gauge 

U{a,t) = 2ap.at , (3.1) 

{Pv = p" the light-cone momentum) the dynamics of the transverse string coordinates 
X°'[a,t) in the plane wave metric 

ds^ = 2dudv + Aab{u)x'^x^du^ + dx^ , (3.2) 
is governed by the quadratic light-cone Hamiltonian 

Hic{t) = -pu = ^^jj- r da (X"2 + x'-2 _ 4a'^plAabi2a'p,t)X''X') . (3.3) 
Note that precisely for the singular plane waves with profile 

Aab{u) = -Ujl6abU~^ (3.4) 

^Such an identification has been used in to define an analogue of temperature for non-equilibrium 
systems. 

^Seejll for a derivation of the Lewis- Riesenfeld theory of invariants of time-dependent harmonic oscil- 
lators |30| from the geometry of plane waves. 
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the dependence of the light-cone Hamiltonian on disappears (up to an overall factor). 
This is a consequence of the scale invariance 

{u,v) ^ {Xu,X-^v) (3.5) 

characterising the corresponding homogeneous plane wave metric |jl|. 

Thus the Fourier modes of the string are harmonic oscillators with frequencies (up to 
a standard overall factor of l/a'p^) 



u:l{t) = ^n^ + ^ . (3.6) 

We will consider the case that all the frequencies Ua are real. This excludes vacuum plane 
waves, but includes e.g. the Penrose limits of FRW metrics Q and a large class of other 
power-law singularities [Q].^ 

The natural mode expansion in terms of solutions to the classical equations of motion 
and the corresponding invariant oscillators (Appendix A) leads to a non-diagonal light- 
cone Hamiltonian operator, namely (cf. |jl2| and ( |A.7] )) 

^ ^ oo d 

HUt) = -f-[Ho{t) + Et^n W(«-n«n + «-n«n) " i^n(i)«n«n " J] 

O- Pv ^ , , 

71=1 a=l 

(3.7) 

where Ho{t) is the Hamiltonian of the zero modes and the time-dependent coefficients are 
quadratic expressions in the Bessel function solutions to the string mode equations. 

As shown in this specific example in [12|, and discussed more generally in Appendix 
A in a quantum mechanical context, the light-cone Hamiltonian operator simplifies signif- 
icantly when written in an explicitly time-dependent but diagonalising basis of oscillators 
a^{t) , d^{t) , in which the coefficients are just the appropriately normalised classical fre- 
quencies (|3.6D rather than the complicated functions il$^(t) appearing in (^?7|), 

1 CO d afj.\ 

a Pt, ^-^ ^-^ n 

n=l a=l 

Here /i(t) is a normal ordering c-function. 



d 



a=l n=l 

which we will discuss in more detail in section 5. The oscillators are normalised (as are the 
a oscillators in ( |3.7D ) to have the non-zero commutation relations 

K(i),a^(t)] =n5n+m<5'^' 

K(t),a^(t)] = nbn^m^'''' . (3.10) 

In particular, in this basis the invariant 

I = Hi,{tQ) (3.11) 

takes a very simple form and its corresponding "thermal" partition function is easy to 

calculate (see section 5). 

*For a dicussion of some aspects of "imaginary" frequencies in plane waves see e.g. |3^ . 
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4. Implications of Scale Invariance 



As we have seen, the string mode frequencies in the scale invariant plane wave are the 
p^-independent but t-dependent 



^^(t) = ^n2 + ^ • (4.1) 

This highlights the special feature of the scale invariant plane waves, namely that for fixed 
t their large n behaviour is exactly that of flat space, 

(4.2) 

The above behaviour of strings in scale invariant homogeneous plane waves should be 
contrasted with the behaviour of the string modes in the plane waves with constant wave 
profile 

Aabiu) = -ulSab ■ (4.3) 

In this case one finds (up to the same overall factor as in (^^)) the now t- independent but 
p^-dependent frequencies [16, ^ 



^niPv) = + a'^PWa ■ (4.4) 

While these also behave as ~ n for large n and fixed p^, the integration over the light- 
cone momentum that arises in the calculation of the partition function implies that for the 
symmetric (i.e. constant Aab) plane waves, or any u-dependence of the frequencies other 
than times a bounded function of n, there can be significant departures from the fiat 
space behaviour even for n — > oo. 

This seems to indicate that the scale invariant singular homogeneous plane waves have 
properties more in common with strings in flat space than either have with the symmetric 
plane waves. Another manifestation of this is the fact that for u > the Wa ^ limit of 
the metric with profile (|3.4| ) is the flat metric, while rescaling Ha in ( f4.3D is an isometry of 
the metric (so that the flat space emerging at /u^ = is a non-Hausdorff limit). 

In particular, this has direct implications for the issue we want to study in this paper, 
namely the Hagedorn behaviour of strings in scale invariant plane waves, which depends 
on the properties of the string modes and the exponential growth of the number of states 
at large n. On the basis of the above reasoning, one might expect this behaviour to be 
identical to that in flat space, and we will confirm this by an explicit calculation in the next 
section. On the other hand, as was already mentioned in the Introduction, the symmetric 
plane waves exhibit a different behaviour. 

5. The "Thermal" Partition Function 

In this section we will present the calculation of the thermal partition function ( ^.6| ) cor- 



responding to the invariant / = Hic{to) (|3.11 ). We will add to the invariant the light-cone 



momentum in such a way that in the flat space limit the density operator reduces to 
the standard expression for the time-like Hamiltonian {pu +Pv)/V^- 
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Thus the density operator of interest is 



= TTr TT > (5-1) 

tre-A(^ + ^^) 

and we will now study the associated thermal partition function 

Zj(/3) =tre"A^^'^ + ^^ . (5.2) 

To evaluate this trace, we need to implement the level matching condition. We first do the 
calculation for a plane wave solution with non-trivial dilaton and no other non-trivial fields 
in the bosonic string theory |jl^ and will discuss later the generalisation to other string 
theories and field configurations. 

In terms of either of the oscillator bases the level matching condition is simply the 
difference between the left and right number operators, 

d oo d oo 

N -N = ^ ^(a" - a" „<) = ^ ^{a^a^ - a'^cQ ■ (5.3) 

a=l n=l a=l n=l 

Thus Zi{(3) is 

Zm- f<lp, /dAe-W..e-^^ + ^"^<'^-*' . (5,4) 



For notational simplicity we will from now on consider the case where all the frequencies 
are equal, tOa = w, but our final result for the critical (Hagedorn) value of f3 turns out to be 
independent of oj and is hence, in particular, also valid when the frequencies are distinct. 
We will also abbreviate LOnito) = oJn- 
In terms of the complex variable r, 

T = ri + iT2 = A + i — j£- , (5.5) 

the trace inside the integral is 

oo oo 

^^^-2ttt2I +2-KiTi{N -N) = | ^ ^ (-27rr2W„ 27rzrin)m|2d 

=1 m 

oo 

E 



n=l m=0 

oo 

-2TTT2iOQd'm — 27rr2/i(to) 



X 2^e 

m=0 

oo 



^ -2TTT2h{tQ) 

-27rT2UJn + 2'KiTin' 



= Doo{ruT2; e2^^2(dAo(g) - h^)) (5 g) 

2to 

where DQfi{T;q) is a generalised (massive) theta function (see Appendix B for the defini- 
tion). 
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The normal ordering term 



oo 



h{to) = diY^ ^jn^ + + it) ' (^-^^ 



n=l 



is divergent. Actually /i(to) has two divergences, the first one arising from a sum over n of n 
- a quadratic divergence, the second one arising from a sum over n of 1/n - a logarithmic di- 
vergence. In the present case (generalising bosonic string theory in Minkowski background) 
this quadratic divergence is cancelled by a generalised zeta function regularisation. The 
details of the resummation are given in Appendix B with the result that {q = uj/2tQ), 



Hto) 



i ^ + g2 (5.8) 



2 



= -iV^l^^-^r(-i) (5.9) 

k=l 

= Ao(g) - ^r(-l) , (5.10) 

where Ki is a modified Bessel function of the second kind. 

For q = (Minkowski space) it is easy to see that we get the usual result of zeta 
function regularisation that /i(to) = —d/12 (see Appendix B). However, for g 7^ the 
term proportional to needs interpreting as it is infinite. The source of this infinity is 
the subleading logarithmic divergence ~ C(l) of /i(to) mentioned above, as can also be 
seen directly by expanding and resumming the series for small q by utilising the binomial 
expansion of the square root (|B.9|) . 

This remaining divergence is more subtle, but before we consider it in more detail we 
first make the following observation. Since, in the present context, q is not related to the 
modular parameter, i.e. independent oipv because of scale invariance, the relevant modular 
transformation rule of DQfi{T;q) is ^ 

^o,o(n, T2; q) = Dofl{-n/\r\^ T2/|r|2; g|r|) . (5.11) 

As a consequence T2q'^ is itself modular invariant, 

r2q' = {r2/\r\'MT\f , (5.12) 

and thus any term in /i(to) that is proportional to q^ makes a contribution to ( |5.6| ) that is 
modular invariant all by itself. 

In principle there can be a logarithmic divergence arising in a two-dimensional non- 
linear sigma model though it generally vanishes on shell or can be removed by a singular 
field redefinition as such models are formally (1-loop) scale invariant. In string theory 
one requires that the quantisation preserves modular invariance and divergences are again 
regularised via singular field redefinitions and appropriate subtractions (for a discussion of 



^This is difFerent from the modular transformation for the BFHP plane wave, for which q depends on 



Pv (4.4) and one has Do,o(T;g) = Do,o{ — i/T;q/\T\) instead, of. the discussion in |l8[ 
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these issues in the present context, see For us at the moment it is important simply 

to note that any subtraction that removes the divergence wih be proportional to and 
thus will not destroy the modular covariance of the state sum. 

Putting everything together, and setting d = 24 for the bosonic string, we see that we 
can now write Zj{(3) as 



Zi{P)= / dne-^^2Do,o{n,T2;^)-'^. (5.13) 

Jo ^2 J~l/2 



The potential divergence in this integral arises from the region where p„ — > oo, corre- 
sponding to T2 — > 0, and we can set ri = to determine the leading behaviour of the 
integral. 

To easily determine the behaviour of the integrand one can use the modular trans- 



formation property (5.11) of Dq.o to deduce that for ri = and in the limit T2 ^ the 



transverse partition function behaves as 



^o,o(ri,r2;:^)-2^-exp[-^24Ao(:^r2)] . (5.14) 

T2 ZtQ 



Combining this with the expansion of Aq (Appendix B), 

one finds that the leading behaviour of the integrand as r2 — > is 

exp[-^(167r2a'-/32)] . 

Thus Zi{j3) diverges for P < Pc, with 

= levr^a'. (5.16) 

This is precisely the value Pc = Ph that corresponds to the Hagedorn temperature for 
strings in flat space ||2^. Note that, in particular, this result is independent of to or, 
more generally, of the frequencies uJa determining the plane wave background. It is also 
manifestly independent of tQ or, in other words, independent of the choice of invariant 
/ = H{tQ) determining the thermal ensemble. 

For the type II superstring theories, we have a similar dilatonic background. The 
bosonic contribution to the "thermal" partition function Zj{P) is as above, with d = 8. 
To include the fermionic contribution we note that in the light cone gauge the fermionic 
Lagrangian is given by 

Lf = ij'd.ip' + i>'d,i,' , (5.17) 

so the fermions do not couple to the parameters of the metric, and therefore its contribution 
to Zj{P) will be the same as in a fiat background. It follows that Zf{P) diverges for P < Pc 
where Pc = Ph = 27r\/ 2a' coincides with the inverse Hagedorn temperature of type II 
strings in flat space p9|. 
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For superstring backgrounds that also have non- vanishing form fields one finds the same 
results. For example, consider a type IIB superstring propagating in an homogeneous plane 
wave metric supported by a Ramond-Ramond 5-form background (the singular analogue 



|] of the BFHP background) 

2 



ds^ = 2dudv — '^x^dv? + dx^ 

Fui2?,A = Fu5e78 = 2— . (5.18) 
u 



As in |16, 17|, the GS light-cone gauge action for a string in this background is 

L = Lb + Lp 

Lb = lid+X'^d^X- - ^X''^) 
2 

Lf = i{e^Td+e^ + 0^'d-9^ - 2je^j-m^). (5.19) 

Notice that once again all pt,-dependence has disappeared. The mode equations for the 
fermions can also be solved in closed form (in terms of Bessel functions). Just as for the 
bosonic modes, their large n behaviour uniformly aproaches that of flat space. Thus the 
resulting critical "inverse temperature" Pc is independent of iv and equal to the type II 
Hagedorn temperature of flat space, in contrast to the result for the BFHP (and other 



symmetric) plane waves [18[- p^ . For other (e.g. NS) backgrounds, the string mode equa- 
tions may be more complicated (and difficult to solve in closed form), but the large n 
behaviour will always be as above, leading to the same conclusions regarding the value of 
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A. Invariant vs. Diagonalising Oscillator Bases 

Consider the time-dependent harmonic oscillator with Hamiltonian 

i7 = -/ + ^x^. (A.1) 

For a time- independent harmonic oscillator one can easily find an oscillator basis which 
gives a diagonal Hamiltonian and such that the oscillators themselves are invariant. For 
a time-dependent oscillator one cannot find a basis for which these two conditions are 
simultaneously satisfied. 

It is easy to find invariant oscillators: one defines a and by the "mode expansion" 

x{t) = az{t) +Q^z{t)* (A.2) 
p{t) = m{az{t) +ah{t)*) (A.3) 
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where z{t) is a complex solution to the equations of motion 

z{t) = -uj{tfz{t) (A.4) 
with the Wronskian of z{t) normalised to 

W{z,z*) = z{t)i*{t) - z{tyi{t) = - . (A.5) 

m 

The Heisenberg operator equations of motion imply that indeed a is invariant, in the 
sense that it satisfies (p.3|), and the Wronskian normalisation condition implies that these 
oscillators satisfy the standard canonical commutation relations 

[a,a^ = 1 . (A.6) 

Note that then also any, say, quadratic function of these oscillators with time-independent 
coefficients is an invariant in the sense of 

In this basis the Hamiltonian takes the general form 

H = ^[(aat + at„)(|i|2 + ^2^,\2^ ^ -2 ^ ^2^2) ^ ^t2(^*2 ^ ^2^*2)] (^^7^ 

It is non-diagonal for io{t) not constant, because 

+ uj{tfz{tf) = 2Lo{t)ij{t)z{tf (A.8) 

at 

is then not zero, so a fortiori z"^ + uP'z^ itself cannot be zero. 

For a;(t) = uj constant, the Hamiltonian is diagonal and explicitly the above mode 
expansion reads 



zmuj 



p{t) = ^^^^-ae-'^^ + aU'^^) . (A.9) 

Since the proof that the Hamiltonian in this basis is diagonal is purely algebraic, i.e. does 
not depend on the t-independence of uj, we are thus led to define, in general, an alternative 
oscillator basis a and by 



= /_?!_(ae-^'^W* + ate^^W*) (A.IO) 
y 2muj{t) 



Pit) = ^y^^(-ae— + ate-W*) . (A.ll) 

This indeed defines the oscilators a and a\ as one can solve for them in terms of x{t) and 
p{t), and this determines the non-trivial time-dependence of a(t) and a^^(t). Nevertheless 
one has the time-independent canonical commutation relations 

[a{t),a\t)] = 1 . (A.12) 
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In terms of these oscillators the Hamiltonian is diagonal, 

H(t)=^!^^aa^ + a^a) (A.13) 

and, in particular, the frequency in the diagonal basis is just the classical frequenncy uj{t). 

Explicitly, the time-dependent SU{\, 1) transformation ( = invariance group of the 
oscillator algebra) from a to a is 

a = f{t)a + g{t)a^ , (A.14) 

where 

»(«) = W'5j^e''^(*"(u.(()3(t)' +;=(<)•), (A.15) 



satisfy the SU{\, 1) condition 



\f{t)\^-\g(t)\^ = -^W = l (A.16) 



as a consequence of the condition on the Wronskian of z{t). 
B. Generalised Theta Functions 



The generalised Theta function is given by |33| 



oo 

Ai,f>2(n,r2; q) = e2^^2Ab,(g) J] (i _ e^vrraV (n + hiY + + 27riri(n + 61) - 2mb2^ ^ 

n=—oo 

(B.l) 



where 

,^cc2e^M;f,(2,,rt, (B.2) 

vr p 
p=i 

and Ki is a modified Bessel function of the second kind. Its modular properties are 

^61,62(^,^2; = Dfe2 7-r^;q\T\) = Db^^bi+b^in + l,r2;g) . (B.3) 

In the detailed calculations of this paper we only need DQfi{Ti,T2; q) and therefore just 
Ao(g), whose relation to the normal ordering c- function /i.(to) ( |5.7| ) we will now explain. 
To that end consider F = 2h{tQ)/d, which we write as 

F=Y, Vn'^ + (B.4) 

where q = uj/2to. In zeta function regularisation the first step we take is to use Poisson 
resummation. Then 

/oo 
^yg-27rz%y^2^. (B.5) 
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Rewriting the square root term as an integral, 



we can then write 



r(3/2 - k) n2fc-i/fc! 

n=l fc=0 ^ ' ^ 

^ + r(3/2-fc)fe! ^ 



r(3/2-fe)fc! 



giving the expansion 



(B.6) 



2g^ i^i(27rfcg) 
fc=i 

= 2Ao(g) - ^r(-l) . (B.7) 

Thus we find 

^ = Ao(,) - ^r(-l) . (B.8) 
Alternatively the binomial expansion of F for small is 

oo 

F = g + 2^ Vn2 + g2 



.+2C(-l)+.-c(l)+2f: ™f-^,^ -^ (B.9) 



(B.IO) 
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